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. ABSTRACT 

<D 

| We present an optimized layout of horn positions for the Degree Angular Scale Interferometer 

(DASI) , which provides good coverage of the u-v plane and 3-fold symmetry. We investigate how 
an optimal subspace filtering analysis of a single pointing of DASI could be used to determine 
the anisotropics in the cosmic microwave background over a range of angular scales near 30'. 
■ We discuss the complementarity between the angular scales probed by DASI and NASA's 

MAP satellite at 30GHz and strategies for imaging the sky at these frequencies. Finally, we 
quantitatively assess how well the angular power spectrum can be recovered by deconvolution of 



^| ' linear combinations of squares of the DASI visibilities. 
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P , 1. Introduction 

h 

The study of the Cosmic Microwave Background (CMB) anisotropy holds the promise of answering many 



of our fundamental questions about the universe and the origin of the large-scale structure (see e.g. Bennett, 
Turner & White 1997). The advent of low- noise, broadband, millimeter- wave amplifiers (Popieszalski 199^ ) 



has made interferometry a particularly attractive technique for detecting and imaging low contrast emission, 
such as anisotropy in the CMB. An interferometer directly measures the Fourier transform of the intensity 
distribution on the sky. By inverting the interferometer output, images of the sky are obtained which include 
angular scales determined by the size and spacing of the individual array elements. 



In an earlier paper (White et al. 1998, hereafter WCDH) we outlined a formalism for interpreting 
CMB anisotropies as measured by interferometers. In this paper we extend this analysis to consider an 
efficient method of analyzing the data that would be obtained from a series of uncorrelated pointings of an 
interferometer. In particular we examine what the upcoming Degree Angular Scale Interferometer^] (DASI, 



Halverson et al. |1998| ) experiment may teach us about cosmology. 



DASI is an interferometer designed to measure anisotropies in the CMB over a large range of scales with 
high sensitivity. The array consists of 13 closely packed elements, each of 20cm diameter, in a configuration 
which fills roughly half of the aperture area with a 3-fold symmetry. Each element of the array is a wide-angle 
corrugated horn with a collimating lens. DASI uses cooled HEMT amplifiers running between 26-36GHz 



1 More information on DASI can be found at 



tittp: //astro .uchicago . edu/dasi . 
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with a noise temperature of < 15K. The signal is filtered into ten 1GHz channels. Details of the layout of 
the DASI horns is given below. 



The outline of this paper is as follows. In §E we adapt the formalism of WCDH to deal with real 



and imaginary parts of the visibilities independently, and show explicitly how the formalism automatically 
imposes the constraint that the sky temperature is real. With the issues defined, we describe the configuration 
of DASI in §|[ In §|[ [| we analyze mock data appropriate to single fields from DASI, showing how to construct 
an optimal basis for likelihood analysis. We specifically address the question of optimal sampling on the sky, 
which was omitted from our previous work. This leads naturally to a discussion of Wiener filtered (Bunn, 



Hoffman & Silk 1996) map making and we reformulate the strategy for imaging the sky in this basis, pointing 



out the complementarity between NASAs MAP satellite?] and DASI at 30GHz. We discuss estimates of the 
angular power spectrum which are easy to implement for uncorrelated pointing of DASI in and show that 
the finite field of view does not hamper our ability to reconstruct the angular power spectrum. Finally in §0 
we discuss multi-frequency observations and present our power spectrum results in the context of "radical 



compression" (Bond, Jaffe & Knox 1998) 



2. Formalism 



The reader is referred to our earlier paper ( WCDH ) for a detailed discussion of how to formulate the 



data-analysis problem with an interferometer, plus references to earlier work. We briefly review the major 
elements here. 

Under the assumption of a narrow frequency band and a distant source, the datum measured by an 
interferometer is proportional to the Fourier Transform of the observed intensity on the sky, i.e. the sky 
intensity multiplied by the instrument beam. We label the "primary" beam of the telescope by A(x), with x 
a 2D vector lying in the plane^ of the sky. Every pair of telescopes in the interferometer measures a visibility 
at a given point in the Fourier Plane, called the u — v plane, 

V(u) cx J dx A(x)AT(x)e 2 " u x (1) 

where AT is the temperature (fluctuation) on the sky and u is the variable conjugate to x, with dimensions 
of inverse angle measured in wavelengths. The (omitted) proportionality constant, dB v /dT where B v is the 
Planck function, converts from temperature to intensity. The spacing of the horns and the position of the 
beam on the sky determine which value of u will be measured by a pair of antennae in any one integration. 
The size of the primary beam determines the amount of sky that is viewed, and hence the size of the "map" , 
while the maximum spacing determines the resolution. 

The 2-point function of the observed visibilities is the convolution of the sky power spectrum, S(u, v), 
with the Fourier Transforms of the primary beams. If our theory is rotationally invariant the power spectrum 



is diagonal, S(u, v) = S(u)S(u — v) and on small scales (WCDH) 



2c , s 1(1 + 1) r 
U S{U) ~ 12^- ° l 



for u £ 10 . (2) 

'—2-KU 



2 http:/ /map. gsfc.nasa.gov/ 

3 Since the field of view of the currently operational or planned instruments is small (< 5°) the sky can be approximated as 



flat with excellent accuracy. 
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where the dimensionless Cg are the usual multipole moments of the CMB anisotropy spectrum and £ ~ O^ 1 
is the multipole index. This approximation works at the few percent level for a standard Cold Dark Matter 
model when u > 10 or £ ;> 60. 

The Fourier Transform of the primary beamj^] is the auto-correlation of the Fourier Transform of the 
point response, g, of the receiver to an electric field, A(u) — g-kg{u) and 

A(x) = J du i(u)e- 2 " u x , (3) 

Due to the finite aperture A has compact support. In order to obtain a simple estimate of our window 
function it is a reasonable first approximation to take A equal to the auto-correlation of a pill-box of radius 
D/2 where D is the diameter of the dish in units of the observing wavelength. Specifically 



~M \ 2A * 

A(u) = 

7T 



arc-cos — 



u u\f D 2 



D D 2 



(4) 



if u < D and zero otherwise. If we require A(0) = 1 then this must integrate to unit area, so A" 1 = ir(D/2) 2 , 
or the area of the dish. We show A(u) in Fig. [|a. For now we shall treat a single frequency. Obviously for a 
fixed physical dish the A will be slightly different for different wavelengths. We return to this complication 
in §@. 

In |WCDH we presented the formalism in terms of complex visibility data. However it is easier practically 



to implement the analysis in terms of the real and imaginary parts of the visibility. We write these as 
Vj = V R + iVj . The cosmological contribution to the real and imaginary components is uncorrelated 
(Vf'Vj) = 0. Assemble V R and V 1 into a vector consisting of first the real and then the imaginary parts - 
the signal correlation matrix of this vector takes block-diagonal form. Further the cosmological contribution 
obeys {V R V R ) — ± \V/Vj). It is straightforward to show that the cosmological contribution is proportional 
to 

~ J rfv S(v)A(u t - v) [A(uj - v) ± A( Uj + v)] (5) 

where the ± refer to the real and imaginary parts respectively. (We have dropped the normalization factor 
which converts temperature to flux.) Note that if u,; = — Uj the visibilities are completely (anti-)correlated, 
as would be expected given that V(u) is the Fourier Transform of a real field: V*(u) = V(— u). The factor 
of i out front reflects the fact that the full variance CX = (VfVj) is the sum of the real and imaginary 
components. 

In the case where all correlated signal is celestial, the correlation function of the noise in each visibility 
is diagonal with 

C»=( 2 J^-) 2 -±-6 t] . (6) 



If the noise in the real and imaginary components is uncorrelated, then each makes up half of this variance. 
Here Ub is Boltzmann's constant, T sys is the system noise temperature, r\A is the aperture efficiency, Ad is 
the physical area of a dish (not to be confused with A(x)), rib is the number of baselines^] corresponding to 



4 Throughout wc will use a tilde to represent the Fourier Transform of a quantity. 
5 The number of baselines formed by n r receivers is nj, = n r (n r — l)/2. 
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a given separation of antennae, A„ is the bandwidth and t a is the observing time. Typical values for DASI 
are T sys = 20K, t]a ~ 0.8, dishes of diameter 20cm, n\, — 3 and A„ = 10GHz (in 10 x 1GHz channels). 

We show in Fig. [l] the positions at which DASI will measure visibilities. With 13 horns there are 78 
baselines. Because of the 3-fold symmetry of the instrument each pointing samples V at 26 different |u|. For 
each "stare" DASI is rotated about an axis perpendicular to the baseplate to fill half of the u — v plane as 
shown in Fig. [|. The other half of the u — v plane is constrained by the symmetry V*(u) = V(— u). 



3. The DASI configuration 

The layout of the horns for DASI is shown in Fig. [l]a. This configuration has a 3-fold symmetry about 
the central horn. The positions of 4 of the non-central horns are arbitrary, up to a global rotation, with 
the remaining configuration being determined by symmetry. The configuration shown represents an optimal 
configuration, within the physical constraints, for the purposes of measuring CMB anisotropy. 

As described in the last section, the distance between each pair of horns represents a baseline at which 
a visibility can be measured. Each visibility probes a range of angular scales centered at £ = 2iru where u is 
the baseline in units of the wavelength. The sensitivity as a function of £ is given by the window function 
(see Eq. ^). Since the width of the window function is determined by the size of the apertures, the optimal 
coverage for the purpose of CMB anisotropy is that configuration which spans the largest range of baselines 
with the most overlap between neigbouring window functions. In one dimension such optimal configurations 
are known as Golomb[] rulers (Dewdney |1985| ), and a well known procedure exists for finding them. 

Finding an optimal solution in two dimensions, within the physical constraints, cannot be done 
analytically. We have optimized the configuration numerically. We have searched the 7 dimensional 
parameter space of horn positions (an x and y position for each of 4 horns, minus one overall rotation) 
for the configuration which minimized the maximum separation between baseline distances, while at the 
same time covering the largest range of angular scales. Trial starting positions were determined by a simple 
Monte Carlo search of the parameter space. From each of these positions a multi-dimensional minimization 
was started. Two additional constraints were imposed upon the allowed solutions: no two horns could come 
closer than 25cm, the physical size of the horn plus surrounding "lip" , and all horns etc had to lie completely 
within the 1.6m diameter base plate (leading to a maximal radius of 63.5cm). Our "optimal" solution is 
given m Table @. The baselines run from 25cm to 120cm with the largest gap between baseline distances 
being 6.4cm. 



4. S/N Eigenmodes or Optimal Subspace Filtering 

A quick glance at Fig. [l]b shows that the signal in most of the visibilities will be highly correlated, i.e. the 
apertures have a large overlap. This is shown quantitatively in Fig. [|a where we plot one row of the signal 
correlation matrix. Given a "trial" theory we can perform a change of basis to remove these correlations (see 
§6.2 of WCDH , or Tegmark et al. 1998 for a review of this method). The input theory can be considered 
as a prior in the context of Bayesian analysis, or could be iterated to match the data if so desired. For our 
purposes all that will matter is that the signal variance is independent of u and decreases with |u|. For 



tittp: / /members. aol.com/golomb20/ 
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concreteness we will take u 2 S(u) ^constant, normalized to the COBE 4-year data (Bennett et al. 1996). 
For this choice the cosmological signal is approximately equal to the noise in each of the highest |u| bins in 
Fig. [|b. 

Let us consider only the real parts of the visibilities for now, the imaginary parts are dealt with in 
an analogous manner. Take these to lie in the upper-left block of the block-diagonal correlation matrix. 
Denote the noise in the real part of each visibility tjj. The eigenvalues of the matrix C RR /<Ti<Tj measure the 
signal-to-noise in independent linear combinations of the visibilities. The independent linear combinations 
can be written v a = J2ii^i/ <7 i) 1 ^'ia where is the ath eigenvector of C RR /ai(Tj. Then {v a Vb) = (A Q + l)5 a b 
where the A a are the signal-to-noise eigenvalues and the 1 represents the noise contribution (which is the 
unit matrix in this basis). 

The number of modes v a with A a > 1 is a quantitative indication of how much signal is being measured 
by the experiment. We show the eigenspectrum of the real parts of the visibilities for 1 day of observation 
in Fig. |b. There is in addition an eigenvalue of the same size for each imaginary component. Due to the 
overlapping windows in the u — v plane, <J 25% of the modes contain good cosmological information. The 
rest are redundant, containing mostly noise. Note however that in 1 day of observation DASI measures 
nearly 200 high signal-to-noise eigenvectors! We have tested that reducing the number of orientations of 
the instrument, i.e. the oversampling in 9 U , does not significantly reduce the number of high signal-to-noise 
eigenvectors until the apertures in the outermost circle are just touching, at which point there are 160 modes 
with S/N > 1. 

The A a are the variance of the v a and can be predicted given a theoretical power spectrum 

A a = \j vdv S(v)W a (v) (7) 

with the window function 

W^(|v|) = V— [d0 v A(u i -v)\A(u J -v)±A(yL j + vj\ (8) 

and as such could form a basis for "radical compression" . We shall return to this issue in WA. 



5. Imaging the Sky 



DASI provides coverage over a significant part of the u — v plane, and as such is able in principle to 
perform high resolution imaging of the sky. The high signal to noise of the DASI instrument and the large 
dynamic range in angular scale however make imaging a computationally challenging task. In this section 
we present two approaches which we have tried with mixed success on simulated data. 



We discussed the Wiener filtering formalism for producing an image from the visibility data in WCDH| . 
Here we note that in the v a basis it is straightforward to construct the Wiener filtered sky map (Bunn, 



Hoffman & Silk 1996; WCDH and references therein). Recall that for Wiener filtering the sky temperature 
is approximated by 



T™(x a ) = C%,W Pj [C v + C 
where C T is the real-space temperature correlation function, 



ijk 



(9) 



(10) 
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we choose cos for the real components and sin for the imaginary components, and C v and C N are the visibility 
signal and noise correlation matrices. We have written the expression in this mixed way to avoid needing 
to invert a matrix of size N pix x N vis where N p - lx is the number of pixels in the map being reconstructed, 
but conceptually the C T W term can be replaced with W~ X C V . Then in the v a basis one simply replaces v a 
with A a (A a + Y)~ x v a which down- weights modes with low signal-to- noise. Specifically 

a a 

where 

M 0a =Y,— A ( x /3) { } [2* * Uj ■ x p ] (12) 

(T i S1H J 

j J 

where we choose cos for the real components and sin for the imaginary components. The ratio A a (A Q + 1) _1 
is plotted in Fig. ||b, where we can see that <J 25% of the modes will contribute significantly to the final 
map. The expected variance in the final map is also easily computed: 

{t wf t w F) = m^-JL-m- 1 (13) 

a a 

which approaches {T 2 ) as A a — * oo. Note that the Wiener filter is not power preserving, so the maps should 
not be used for power spectrum estimation. 

By reducing the total number of modes that need to be kept in the summation the v a basis can speed 
up calculation of the Wiener filtered map. This is an advantage, but even with this speed up we found that 
producing a high resolution map of even a single field on the sky is very computationally expensive due to 
the large number of matrix multiplications involved in computing T . An alternative to Wiener filtering, 
which is very similar in the high signal-to-noise regime in which we are working, is the minimum variance 
estimator for T(x). The number of operations required to produce the minimum variance and Wiener filtered 
maps are comparable, and both tend to be very slow. We are currently investigating faster approximate 
methods of image making. 

The formalism above can in principle be used for making maps at each frequency, or one can make a 
map which combines the frequencies in such a way as to isolate the CMB (or foreground) signal. While 
the formalism looks more complex, it is in fact easy to implement computationally. The techniques are well 
known (see we imagine that at each point in the u — v plane our visibilities form a vector V whose 
components are the different frequencies. We can expand this vector in terms of the physical components 
we wish to consider, e.g., 

V = O C MBV CMB +0ffV ff (14) 
In the absence of noise it is easy to solve for #cmb as 

0cmb = with V± ■ Vff = (15) 

V± ■ VqmB 

If we write #cmb = ^2 a c a V a then we should replace A by J2 a c a,A a throughout. In the presence of noise we 
may define c a by a minimum variance estimator as above: 

c« = £fen CMB ^V| ( 16 ) 
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where N a b is the channel noise matrix and capital Roman letters indicate the component being considered 
(see ©. 

We end our discussion of imaging with an observation about the complementarity between the MAP 
satellite and DASI. In making maps of the sky with DASI, the largest source of error is the missing long 
wavelength modes which are filtered out by the DASI primary beam. The long wavelength power can be 



included in the map by simultaneously fitting another data set (e.g. WCDH ). In this particular case the data 
from the MAP satellite at the same frequency is the obvious choice since MAP will have full sky coverage. 
We show the window functions for MAP at the relevant frequencies, along with the envelope for DASI in 
Fig. |. 



6. Power Spectrum Estimation 

While the v a are the natural basis from the point of view of signal-to-noise and Wiener filtering, and 
can dramatically improve likelihood function evaluation, they are not necessarily the quantities of greatest 



physical interest for power spectrum estimation. In WCDH we discussed estimating a series of bandpowers 



using the quadratic estimator of (Bond, Jaffe & Knox 1998] , Tegmark 1997). Here we present a simpler 
approach more appropriate to unmosaiced fields. In this presentation we shall assume we are dealing with 
single frequencies; the case of multiple frequencies is dealt with in 

In CMB anisotropy observations with single dish experiments much has been written about optimal 
ways to estimate the power spectrum. The principle reason is that care must be taken in weighting the data 
to ensure that no sharp cut-offs in real space are introduced. These lead to ringing in Fourier space and 



delocalize the window function (e.g. discussion in Tegmark 1996 ). For the interferometer no such problem 
arises - each visibility samples a compact region in u. For each Vi, the square is a noise biased estimate of 
the power spectrum convolved with a window function. If we define Si = 2 (V^ — No) with N# the noise 
variance in (the component of the) visibility i, then Sj is an unbiased estimator for C^; 

udu S{u)Wu{u) (17) 



where Wu(u) is the window function in analogy to Eq. (g). We have included the factor of 2 in the definition 
of Si to account for the fact that each component Vi of the visibility contributes half of the total variance of 
|Vi| 2 . Take a weighted average of the s»: Sa = E A iSi- The simplest weighting is to sum all of the Si with 
~ u and we shall use that below. Under the assumption that the visibilities are Gaussian, the error 
matrix for the estimates Sa is 

(SS A SS b ) =2Y J E A i{Cl j +C!>fE jB . (18) 

ij 

Note that for N uncorrelated visibilities (real and imaginary parts) contributing to Sa this gives SSa/Sa = 
7V _1 / 2 (1 + noise/signal) as expected. 

For DASI in the configuration shown in Fig. [l]b one can construct 26 different estimates S A , which will 
however be quite correlated. We show the 13 x 13 correlation matrix for every second estimate in Table || 
along with the expected error on each determination (the diagonal elements). As independent pointings are 
included in the analysis the error bar on each S A decreases as A -1 / 2 , but the correlations remain the same. 



Without mosaicing (WCDH) the resolution in £ is restricted to ~ 2ttD, thus the individual determinations 



are required to be highly correlated. The error on the highest i bins is still dominated by the small number of 
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independent samples, in this case 16. Increasing the oversampling in the angular direction (and the observing 
time) can reduce this error to about 23%, hardly worth the extra time compared to including a different 
pointing. 

The Sa are estimates of the power spectrum convolved with the window function Wa(u). Given a theory 
it is straightforward to compare to the data once the window functions are known, and we discuss this in §|?]. 
However one could ask whether (or how much) cosmological information has been lost due to the convolution 
or if the information is still present in the correlated visibilities. One way of answering this question is to 
attempt to perform an approximate (theory independent) deconvolution. 

The simplest ("direct") method is to assume that S(u) ^constant through the window, so that S(u) » 
Sa/ Judu Wa{u) with Wa{u) — Yli^AiWu(u). (Alternatively one could assume that u 2 S(u) ^constant, 
which leads to a similar expression.) For the window function of Eq. (jj) we have J d 2 uA 2 (u) ~ 0.585Z?~ 2 . 
As we show in Fig. pi this technique works surprisingly well for DASI, even though the FWHM of each 
window is A£ ~ 100, over which scales we can expect power spectra to change significantly. For a CDM 
power spectrum for example, the approximation above induces a systematic 10% error (at worst) due to the 
window function "washing out" the peak structure. 

To avoid this problem, we can attempt to perform the deconvolution by an iterative procedure. Recall 
that we are trying to constrain the power spectrum, while our measurements are the spectrum convolved 



with a positive semi-definite kernel - the window function. In this situation Lucy's method (Lucy 1974) can 



be used. We have implemented Lucy's algorithm, following Baugh & Efstathiou (1993), including a coupling 
between different bins and a regularization of the iteration. As they found, the final result is not sensitive 
to the mechanism chosen. 



To briefly re-cap the method: we think about the deconvolution problem, following Lucy (1974), first 
in terms of probability distributions. If we denote by p{x) the probability of measuring a quantity x and 
p(y\x) the conditional probability of measuring y given that x is true then 

P(y) = J p(y\x)p(x)dx (19) 

which is a convolution integral. We wish to estimate p(x) given observations p ohs (y). We start the rth 
iteration with an estimate p r {x) of p{x) and predict p r {y) using Eq. (jl9|) assuming p{y\x) is known. Writing 



the inverse of Eq. (19), for p(x), using the observed p (y) and rewriting p{x\y)p{y) — p{y\x)p{x) leads us 



to the iterative expression for p r+1 {x): 

P + (x) = P (x) — r , i ' (20) 

J P{y\x)dy 

where the denominator is unity. The iterative method we use takes this expression over with the replacements 
p(y) — > Sa, p(x) — > fiu) = u 2 S(u) and p{y\x) — > u Wa(u). Approximating the integrals as sums equally 
spaced in u we have the iterated pair of equations: 

Sa = ^fK)^ 1 ^^) Ah (21) 

r+ , K) . ; w s(|wj (22) 

£ 2 Ua W A {u a ) 

To make the iteration converge more stably we in fact replace only a fraction e of f r with f r+1 on each step, 
and average the f r {u a ) using a 2nd order Savitsky-Golay filter of length (2, 2). The final result is insensitive 
to the details of this procedure, the number of bins chosen for u a etc. 
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We show in Fig. |5| results of a deconvolution assuming that 16 independent patches of sky were observed, 
each for 1 day. We can see that the deconvolution works well, with good resolution in £ and no bias over 
the scales where DASI has sensitivity. Thus we expect that no important cosmological information has been 
lost by the convolution procedure. The points are highly correlated, but we do not show the correlations on 
the Figure (all correlations were included in the analysis). The error bars are so large because they are the 
computed from the variance allowing all other bins to vary freely. 

Though the example here was for a single frequency channel, the generalization to multiple frequencies 
is straightforward and can in principle allow even finer sub-band resolution. 

Finally we remark that obviously a statistical comparison of a given model to the data should be done 
with the Sa- The correlations of the Sa can be computed for any given theory (e.g. Table ||). This allows a 
full likelihood analysis to be done, and is the route which should be taken when comparing DASI data to a 
specific theory. This is what we shall discuss now. 



7. Radical Compression 

To constrain theories using interferometer data is much easier than in the case of single dish data, thus 
many of the powerful techniques developed for the latter case are not needed. The primary reason for this 
is that interferometers work directly in the £-space of theories. It is straightforward to perform "radical 



compression" (Bond, Jaffe & Knox 1998) of an interferometer data set and quote a set of bandpowers along 
with their full noise covariance matrix and window functions. We shall develop this idea briefly in this 
section. 

We now reintroduce the multi-frequency nature of the data set that has been suppressed during most 
of this paper. For DASI we can work with 10x1 GHz channels which we shall label with a greek subscript. 



Following Dodelson (1997; see also White 1998) we imagine that our visibility signal V a at each point in 
the u-v plane is a sum of contributions with different frequency dependences: V a = Xm^^q 4 - Let A = 
be the CMB contribution, whose frequency dependence will be V° — (1, 1, • • •) for observations at DASI 
frequencies. At each frequency, a, and visibility position, i, the signal is the convolution of the sky with an 
aperture Af which will be of the form of Eq. (Q) with the central u and D varying by the inverse of the 
observing wavelength. If each visibility has noise we can estimate the CMB component 9° by minimizing 

X 2 = E E (^° bs - A V a A ) N£ (n° bs - B V b B ) (23) 

ab AB 

where we have suppressed the visibility index i. Solving d\ 2 j d9 A = amounts to taking a linear combination 
of the frequency channels 9® — c aKi 0bs with 



^=e(e^x- c V) e^ 1 

A V be Id 



(24) 



While formidable this expression reduces to the well known least-squares weighting in the limit N a b = o-^5 a b: 

Ab ^ a b J a 

Now we simply replace Vi with 9® and Ai with ^ a c a Ai a to generalize Eq. (p"7|). The generalized Sa 
can be easily calculated from the data, and the expectation values and distribution can be calculated for any 
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thcory once the window function and noise properties are given. Supplying this set would be an ideal way 
to release the DASI data for the purposes of model fitting. 



8. Conclusions 



We have shown how one can implement the formalism of WCDH in the case of a single stare of the DASI 
instrument. We have discussed making maps, filtering the data and reconstructing the anisotropy power 
spectrum. Our results suggest that with ~ 1 month of data, DASI could provide significant constraints on 
theories of structure formation. 



We would like to acknowledge useful conversations with Alex Szalay and thank Tim Pearson for 
comments on a draft of this work. 
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Fig. 1. — (a) The positions of the dishes for DASI. Each dish has an effective optical diameter D = 20A with 
A = 1cm (solid lines). The physical size of the dishes is 25cm (dotted lines) and the base plate for DASI is 
1.6m in diameter (dashed line), (b) The positions at which visibilities will be measured in 1 "stare" of DASI. 
The circles show apertures of radius 10 (roughly the half power point of the beam profile). Each visibility 
consists of an independent real and imaginary part. Solid lines indicate where data will be taken, the dotted 
lines indicate points constrained by symmetry of the transform. 




Fig. 2.— (a) The aperture function A described in the text (Eq. |). For DASI, D = 20. (b) The window 
functions for the 26 baselines of DASI. 
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Fig. 3. — (a) One row of the signal correlation matrix, plotted in the u — v plane, (b) the signal-to-noise 
eigenvalues A a for the real part of the visibilities. For every real A a there is an independent imaginary mode 
with the same signal-to-noise. All calculations assume a scale-invariant COBE normalized spectrum. 




Fig. 4. — The window functions for MAP and DASI, showing the complementarity in angular scale of the 
two instruments at low frequency. The MAP data will provide large-scale information over the whole sky 
which can be used when making maps with DASI, which is sensitive to the small angular scales which MAP 
cannot resolve. 



- 13 - 



Direct 




200 400 600 

i!=2ttu 



Fig. 5. — Estimates of the power spectrum obtained from simulated data from 16 independent DASI fields, 
each observed for 1 day. (Upper) The open circles indicate estimates obtained by assuming S(u) ^constant 
through the window function. (Lower) The solid squares are from Lucy's method. The points are highly 
correlated (see text). 



x (cm) 


V (cm) 


+00.00 


+00.00 
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+00.00 
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+34.25 


+36.98 


+08.04 
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+29.18 


-25.45 


+28.01 
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-30.38 


-52.62 


+02.93 


-63.43 


-11.53 


-36.05 


+43.79 


-36.57 



Table 1: The positions of the centers of the horns on the DASI base plate for the optimal configuration 
described in §||. 
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a 


0.49 


0.40 


0.37 


0.33 


0.31 


0.29 


0.28 


0.27 


0.26 


0.26 


0.26 


0.26 


0.26 



Table 2: The Sa correlation matrix, Cab/°~a&b with a a = C/ A , for every second estimate of the power 
spectrum from a single pointing of DASI. Correlations are listed as percentages. The final row shows the 
relative error on the power spectrum in the bin. Since the error is dominated by sample variance, for N such 
pointings a is reduced by TV 1 / 2 . All calculations assume a scale- invariant COBE normalized spectrum. 



